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ABSTRACT 

The age calibration of the Washington 5T\ index is mainly used to estimate ages 
of star clusters older than 1 Gyr, no age-metallicity degeneracy effect is considered. 
We have profusely exploited synthetic T) versus C — T) colour magnitude diagrams 
aiming at exploring the intrinsic behaviour of the 8T\ index. The analysis shows that 
STi varies with age and metal content as well. In general, the dependence on age 
weakens for ages greater than ~ 6 Gyr, and results even less sensitive to age as the 
metallicity decreases. For ages younger than ~ 5 Gyr ST\ shows a strong correlation 
with both age and metallicity. The SC index -defined as 8T\ for the C passband- is also 
a combined measurement of age and metallicity. We introduce a new age-metallicity 
diagnostic diagram, 5T\ versus SC - 5T \, which has shown the ability of unambiguously 
providing age and metallicity estimates, simultaneously. The new procedure allows to 
derive ages from 1 up to 13 Gyr and metallicities [Fe/H] from -2.0 up to +0.5 dex, 
and is independent of the cluster reddening and distance modulus. It does solve the 
constraints found in the STi index and surpasses the performance of the standard 
giant branch metallicity method. All these features make the diagnostic diagram a 
powerful tool for estimating accurate ages as well as metallicities. 
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1 INTRODUCTION 


|Geisl er et al. (1997, hereafter G97) provided a calibration for 
the magnitude difference (5) between the giant branch clump 
in intermediate-age star clusters (the horizontal branch in 
old clusters) and the main sequence turnoff as a function 
of the cluster age for the Washington photometric system. 
Particularly, they used the Ti versus C — Ti colour mag¬ 
nitude diagra m (C MP) ( effect ive wavelengths: C ~ 3900A, 


Ti ~ 6300A; Canterna (1976)) of six known star clusters 


(ages > 1 Gyr) in order to measure <5Ti and fitted those val¬ 
ues with the clusters’ ages (see their Fig. 5 and eq. 4). STi 
emulated the former SV index defined by j Phelps, Janes fc| 


Montgomery (1994). By using such a calibration G97 firstly 


searched for old star clusters in the Large Magellanic Cloud; 
the 25 candidates analysed resulted to be of intermediate age 
(1-3 Gyr). Since then, the 8T\ index was employed succes¬ 
sively to estimate the ages of star clusters in the Milky Way 


(e.g. Piatti, Claria & Ahumada 2004), in both Magellanic 
Clouds (e.g. Bica et al.||1998 Piatti et al.||200l I and to de¬ 


rive ages for the so-called representative stellar population 
of galactic star fields (e.g. Piatti, Geisler & Mateluna|20121, 
among others. 

The age calibration of the 8T\ index is based on CTi 
photometry of NGC2213 and ESO 121-SC03 in the Large 
Magellanic Cloud and the Galactic open clusters NGC 1245, 
Tombaugh 2, M67 and NGC 6791. At the very least the cal¬ 
ibration needs revision, since updated ages and metallici¬ 
ties for the calibration clusters are available. For instance, 


Anthony-Twarog, Twarog & Mayer (2007) derived an age 
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of (7.0 ± 1.0) Gyr for NGC 6791 QFe/H]=(+0.42±0.05) dex 
Heiter et al. 2014), whereas a value of 10 Gyr was used by 
G97. Except for ESO 121-SC03 ([Fe/H]=(-0.93±0.20) dex 
Olszewski et al. 1991), all of these calibration clusters are 
considerably more metal-rich than [Fe/H] ~ -0.5 dex. In ad¬ 
dition, the solar abundance cluster M 67 plays a pivotal role 
in the calibration as it is the only cluster in the 2-9 Gyr 
range. Secondly, thinking in the impact that the adopted 
cluster fundamental parameters can have on the former age 
calibration of the 5T\ index, it seems necessary to enlarge the 
cluster sample in order to support a robust relationship and 
to allow an analysis about the possible metallicity sensitivity 
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and the effect any such sensitivity might have on the derived 
ages. Nowadays, however, a representative sample of clusters 
with age and metallicity values distributed well along the 
known cluster age/metallicity regime and with Washington 
CT\ photometry is not available. Given the usefulness that 
5Ti has shown in the literature as a cluster age indicator, 
the above constraints seriously blur its full scope. 

Fortunately, synthetic cluster CMDs are powerful tools 
to probe the genuine performance of STi as an age indicator 
and to disentangle any metallicity dependence. This is be¬ 
cause they can accurately reproduce CMDs of clusters of any 
age and metallicity, bearing in mind the fraction of binaries, 
the cluster initial mass function, the cluster mass, as well 
as photometric uncertainties and completeness effects, etc 


(Popescu, Hanson & Elmcgrccn 2012 Casagrande & Van- 
denB erg|2014[ ). In this sense, synthetic cluster CMDs have 
the advantage of representing those of calibration clusters 
with ages and metallicities derived on a homogeneous scale. 
Indeed, we have taken advantage of the recently developed 
Automated Stellar Cluster Analysis package (ASteCJpl P er-| 
ren, Vazquez & Piatti|2015 \, as we describe in Sect. 2 of this 


work, to profusely exploit a large number of synthetic cluster 
CMDs. From the outcomes of this comprehensive analysis, 
we present in Sect. 3 a new age-metallicity diagnostic dia¬ 
gram that involves S values for the C and Ti bandpasses. The 
proposed technique allows to estimate ages and metallicities 
for clusters older than 1 Gyr, independent of their redden¬ 
ing and distance moduli. Sect. 4 analyses the performance 
of the new procedure in the light of published accurate age 
and metallicity values as well as of different age and metal¬ 
licity calibration for the Washington system. Finally, Sect. 
5 summarises our main results. 


2 SYNTHETIC COLOUR-MAGNITUDE 
DIAGRAMS 

Synthetic stellar populations have usually been generated to 
study the astrophysical properties of different stellar popula¬ 
tions through the comparison of their respective CMDs (e.g. 
Romeo et al.||1989[ |Hernandez, Valls-Gabaud fe Gilmore| 


1999 Monteiro, Dias fe Caetano|2010|. The recovery of the 


star formation history of galaxies ( Rubele et al.||20l2 i, the 

study of the extended main sequence turnoff phenomenon in 


star clusters (Correnti et al. 20141, and the estimation of star 
cluster parameters ( Donati et al.||2015 l have been topics of 
astrophysical interest approached by synthetic CMD anal¬ 
yses. Likewise, synthetic CMDs offer the possibility to ex¬ 
plore in detail the behaviour of different photometric indices 
in terms of astrophysical quantities (Catela n et al.| 2014l, to 
predict the range of distinct fundamental properties in star 
clusters ( Popescu fe Hanson|2014 l as well as to assess in ad¬ 
vance the performance of astronomical instruments (Kerber 


et al. 2009). 


In order to revise the STi age calibration and its possi¬ 
ble dependence with the metal content, we have employed 
the ASteCA suit of functions to generate synthetic CMDs 
of star clusters covering ages from 1.0 up to 12.6 Gyr and 
metallicities in the range [Fe/H] = (-2.0 - +0.5) dex. ASteCA 


http://asteca.github.io/ 


was designed as a new set of open source tools for an objec¬ 
tive and automatic analysis of large cluster data sets. The 
code includes functions to perform cluster structure analysis, 
luminosity function curves, integrated colour estimates sta¬ 
tistically cleaned from field star contamination, a Bayesian 
membership assignment algorithm, and a synthetic cluster- 
based best isochrone matching method to simultaneously 
estimate clusters’ properties (age, metallicity, distance, red¬ 


dening, mass and binarity). Perren, Vazquez & Piatti (2015) 
showed that it does not introduce any biases or new correla¬ 
tions between the various derived cluster parameter values. 

The steps by which a synthetic cluster for a given set 
of age, [Fe/H], distance modulus m — M, and reddening 
E(B — V ) values is generated by ASteCA is as follows: i) a 
theoretical isochrone is picked up, densely interpolated to 
contain a thousand points throughout its entire length, in¬ 
cluding the most evolved stellar phases, ii) The isochrone is 
shifted in colour and magnitude according to the E(B — V) 
and m — M values to emulate the effects these extrinsic 
parameters have over the isochrone in the CMD. At this 
stage the synthetic cluster can be objectively identified as 
a unique point in the 4-dimensional space of parameters 
( E(B — V ), m — M , age and metallicity). iii) The isochrone 
is trimmed down to a certain faintest magnitude according 
to the limiting magnitude thought to be reached, iv) An 
initial mass function (IMF) is sampled in the mass range 
[~0.01 —100] M 0 up to a total mass value M to tai provided 
via an input data hie that ensures the evolved CMD regions 
result properly populated. Currently, ASteCA lets the user 
choose between three IMFs (Kroupa, Tout fe Gilmore|1993 


Chabrier 2001 Kroupa 2002) but there is no limit in the 


number of distinct IMFs that can be added. The distribu¬ 
tion of masses is then used to obtain a properly populated 
synthetic cluster by keeping one star in the interpolated 
isochrone for each mass value in the distribution, v) A ran¬ 
dom fraction of stars are assumed to be binaries, which is set 
by default to 50% (von Hi ppel|2005 1, with secondary masses 
drawn from a uniform distribution between the mass of the 
primary star and a fraction of it given by a mass ratio param¬ 
eter set to 0.7. Both quantities can be modified through the 
input data hie. vi) An appropriate magnitude completeness 
and an exponential photometric error functions are finally 
applied to the synthetic cluster. 

As for our purposes, we used the theoretical isochrones 
computed by jBressan et al.| ( |2012| ) using extensive tabula¬ 
tions of bolometric corrections with uncertainties ~ 0.001 


mag for the C and Tj hlters, the IMF of Chabrier (2001), 


and a cluster mass as a function of the cluster’s age given 
by the expression log(M tota i M@) = 1.8xlog(age yU 1 ) 


12.8 ( Baumgardt et ak||2013 |de Gr ijs, Goodwin & Anders 

[2013; jPiatt i 2014J) in order to keep not only the main se¬ 

quence turnoff (MSTO) but also the red clump (RC) sim¬ 
ilarly populated as those of known low mass clusters. We 
considered no binarity effect, since it has been shown that 
binary stars can broaden the MSTO region (Li et ak||2012 


Li, de Grijs & Deng 2013 Jiang, Han fe Li|2014 |, in contrast 
with our aim of sizing up the sensitivity of the MSTO to age 
and metallicity. Fig. 1 shows some examples of the result¬ 
ing synthetic CMDs. Note that the theoretical isochrones of 
Bressan et al.| ([2012]) are in a very satisfactory agreement 
with another known set of isochrones for the Washington 
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system computed by Lejeune & Schaerer (20011 (e.g. Piatti 
et al. 2003 and references therein). 


3 THE <m INDEX 

We measured the difference in magnitude between the 
MSTO and the RC in the generated synthetic CMDs fol¬ 
lowing the precepts of Phelps, Janes & Montgomery! (1994, 
see their Fig. 1), as G97 also applied in their definition of 
5Ti. The synthetic cluster CMDs allowed us to estimate the 
errors involved in measuring 5T\. We found that an uncer¬ 
tainty between 0.05-0.10 mag typically affects the measure¬ 
ments of 5Ti, and can reach a few more hundredth of mag¬ 
nitude for some of the youngest clusters in our synthetic 
sample because of the shape of their MSTOs. Fig. 2 (upper 
left-hand panel) shows the resulting STi values as a func¬ 
tion of age for six different metallicity levels (coloured line 
scale). As can be seen, the relationship of ST] with age varies 
with the metal content as well. In general, the trend of 8T\ 
with the age weakens (it changes to a smaller slope) for ages 
larger than ~ 6 Gyr, and results even less sensitive to age as 
the stellar population becomes more metal-poor ([Fe/H] < 
-1.0 dex). This means that 5T\ does not result a good indica¬ 
tor of age independent of metallicity for old and metal-poor 
clusters. For ages younger than ~ 5 Gyr 5T\ shows a strong 
correlation with both age and metallicity. 

We included in Fig. 2 (upper left-hand panel) the rela¬ 
tionship obtained by G97 represented by a solid black line. 
For ages younger than 3 Gyr they used three calibration 
clusters, two of them (NGC 2243, Tombaugh 2) have an av¬ 
erage spectroscopic metallicity of [Fe/H] = (-0.4 ± 0.10) 
dex. Their 5T\ versus age relationship roughly matches the 
synthetic one for [Fe/H] = -0.5 dex. On the other hand, the 
older end of the G97’s curve overlaps those of metal-poor 
clusters ([Fe/H] < -1.5 dex). The upper right-hand panel of 
Fig. 2 illustrates clearer this metallicity dependence of 5T\ 
for twelve different age levels. Indeed, for 5T\ = 3.2 mag an 
age of 10.3 Gyr comes out from the G97’s calibration as well 
as from synthetic CMDs with [Fe/H] < -1.5 dex. However, 
an age ~ 2.5 Gyr smaller results from the synthetic CMDs 
if a [Fe/H] > -1.0 dex is adopted. Such a difference between 
mean age values results statistically significant even if we 
consider an uncertainty of 0.15 mag in the STi values. Sim¬ 
ilarly, for STi = 1.5 mag an age of 1.8 Gyr is derived from 
the G97’s equation, or indistinctly from synthetic CMDs 
with [Fe/H] « -0.6 dex. Once again, if [Fe/H] values of -0.1 
dex and -1.1 dex were instead taken into account, the ages 
from synthetic CMDs would result 0.5 Gyr greater and 0.6 
Gyr smaller than the G97’s age, respectively; errors of 0.2 
Gyr due to uncertainties in STi were estimated. 

The above examples not only illustrate that the STi 
age index is sensitive to metallicity, but also that such a de¬ 
pendence is a complex function of both age and metallicity. 
Moreover, the age range for any particular 5T\ value ob¬ 
tained from the synthetic curves (see upper left-hand panel 
of Fig. 2) results many times larger than the age errors de¬ 
rived from eq. 4 of G97. Fig. 3 illustrates this trend for three 
different ST\ error levels. This result points out to the need 
of a new calibration for the STi age index which prevents us 
against any age/metallicity degeneracy. 

We seek for any possible straightforward relationship 


between the ST\ values derived from synthetic cluster 
CMDs, the age and the metallicity. Bearing in mind the kind 
of arithmetic expressions employed in previous 8 age calibra¬ 


tions (Phelps, Janes & Montgomery 1994; Carraro & Chiosi 


1994 Geisler et al.|1997 Salaris, Weiss fe Percival|2004 1, we 

tried different possibilities which included linear, quadratic 
and terms of higher degree in any of the three quantities, 
mixed terms, logarithmic functions, etc. Note that we used a 
larger number of points than any previous 8 age calibration, 
uniformly distributed throughout the whole age/metallicity 
range and without any constraint from non-homogeneity in 
the age/metallicity values. Unfortunately, we did not attain 
any satisfactory fit, which confirms the complex interdepen¬ 
dence of the three parameters. 

In order to complement this analysis and by taking ad¬ 
vantage of the availability of magnitudes in the C passband 
(= Ti + (C — Ti)), Fig. 2 (bottom panels) also depicts the 
relationship for SC - defined as 5T\ but for the C filter - with 
the age and the metallicity. In order to build that relation¬ 
ship, we first produced synthetic C versus C — Ti CMDs for 
the same synthetic clusters used to construct the top panels 
of Fig. 2. Then, we measured the magnitudes at the MSTO 
and the RC and computed their difference (SC). The result¬ 
ing curves show that SC expands over a dynamical range 
of ~ 4 mag, similar to that of 8T\. However, although the 
overall appearance of SC and 8T\ curves look alike, the for¬ 
mer show a less complex trend with age and metallicity. 
Particularly, they do not account for: i) the lack of metal¬ 
licity sensitivity ([Fe/H] > -1.0 dex) for stellar populations 
older than ~ 5 Gyr and, ii) the slope change for clusters 
more metal-poor than [Fe/H] < -1.0 dex. Nevertheless, Fig. 
2 clearly suggests that SC is also a combined measurement 
of age and metallicity. 


3.1 A new age-metallicity diagnostic diagram 


At this point, we decided to introduce a new diagnostic di¬ 
agram with the ability of unambiguously providing age and 
metallicity estimates within certain Washington CTi photo¬ 
metric error limits. We found that the plane 8T\ versus SC - 
STi resulted to be the one which best distinguishes changes 
in age and metallicity throughout the whole two-dimensional 
space. In Fig. 4 we have traced iso-age lines and marked 
iso-abundance positions using colour-coded lines and filled 
circles, respectively. Errobars for typical uncertainties in 
STi and SC - ST\ are also indicated. Note that the iso¬ 
age lines and iso-metallicity positions rely on the theoret¬ 
ical isochrones used to build the synthetic CMDs, while the 
quoted errors refer to the uncertainty in measuring SC and 
5T\ in observed cluster CMDs. In this sense, the capacity 
of resolving ages and metallicity varies with the position 
in that plane. Table 1 lists the age/metallicity errors associ¬ 
ated to typical a(ST\) and a [SC - 8T\) uncertainties, derived 
by interpolation of Fig. 4. Particularly, we have highlighted 
with boldface characters those errors for (age, [Fe/H]) pairs 
that are within the ranges of the age-metallicity relation¬ 
ships (AMRs) of the Milky Way and of the Small/Large 
Magellanic Clouds (|Beasley et al.||2004| |Sabbi et al.||2009| 
Piatti fe Geisler|2013 Since star clusters can have ages and 
metallicities directly linked to the chemical evolution history 
(AMR) of their host galaxies, the highlighted age/metallicity 
ranges in Table 1 result more usable from an astrophysical 
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point of view. We thickened the line sections in the 5T\ ver¬ 
sus age diagram (see Fig. 5) corresponding to age/metallicity 
values with uncertainties in Table 1 in boldface. Thus, read¬ 
ers can consistently compare the G97 locus with the present 
theoretically-driven calibration according to known AMRs. 
Fig. 4 is aimed at entering with ( SC-STi , 5T\) values and 
obtaining by interpolation age and metallicity estimates. 

We used the new diagnostic diagram for clusters 
with high-quality CTi photometry, particularly with well- 
identified MSTOs and RCs, and with accurate ages and 
metallicities. We searched for metallicities obtained from 
high-dispersion spectroscopy, although in some cases we re¬ 
laxed this requirement down to medium-dispersion spec¬ 
troscopy or even to reliable photometric metal abundances. 
As for the cluster ages, we took advantages of those val¬ 
ues derived from isochrones fitting to deep cluster CMDs. 
We excluded any previous age/metallicity estimates coming 
from Washington photometry. 

From barely 340 clusters with published CTi data, we 
found nearly 150 with recognisable MSTOs and RCs. In gen¬ 
eral, the extracted cluster CMDs present signatures of con¬ 
tamination from the composite star field population along 
the line of sight. Such a field contamination does have a par¬ 
ticular pattern given by the luminosity function, the colour 
distribution and the stellar density towards the cluster re¬ 
gion. For these reasons, we first built CMDs representing the 
field along the line of sight towards the individual clusters, 
which we then used to clean the cluster CMDs with the aim 
of tracing fiducial cluster sequences as accurate as possible. 
We employed the cleaning procedure developed by |Piatti fe| 
I S ic~| ( |2012| see their Fig. 12). The method compares the 
extracted cluster CMD to four distinct CMDs composed of 
stars located reasonably far from the object, but not too far 
so as to risk losing the local field-star signature in terms of 
stellar density, luminosity function and/or colour distribu¬ 
tion. Each field region covers an equal area as that of the 
cluster and the four field areas are placed to the north, east, 
south and west from the cluster. 

The comparison of the cluster/field CMDs is performed 
using boxes of different sizes distributed in the same manner 
throughout both CMDs, thus leading to a more meaningful 
comparison of the numbers of stars in different CMD re¬ 
gions than using boxes fixed in size and position. The latter 
is not universally efficient, since some parts of the CMD are 
more densely populated than others. For instance, to deal 
with stochastic effects at relatively bright magnitudes (e.g., 
fluctuations in the numbers of bright stars), larger boxes are 
required, while populous CMD regions can be characterized 
using smaller boxes. Since the procedure is executed for each 
of the four field-CMD box samples, it assigns a membership 
probability (P) to each star in the cluster CMD. This is 
done by counting the number of times a star remained un¬ 
subtracted in the four cleaning runs and by subsequently 
dividing this number by four. For our purposes, we used 
stars that are predominantly found in the cleaned cluster 
CMDs (P > 75%). 

We used the cleaned cluster CMDs to measure C and Ti 
magnitudes at the MSTO and RC, then computed SC and 
STi and entered into the age-metallicity diagnostic diagram 
to estimate cluster ages and metallicities. The resulting val¬ 
ues with their uncertainties for clusters that fulfilled the re¬ 
quirements mentioned above (e.g. available age/metallicity 


values from independent Washington techniques) are listed 
in the last columns of Table 2, in which we also included 
age/metallicity values adopted by thoroughly searching the 
literature for comparison purposes. 


4 ANALYSIS 


From the comparison between cluster ages taken from the 
literature and those estimated above (see Fig. 6) we ob¬ 
tained a mean age dispersion of A(log(age yr“ ) = 0.07 ± 
0.02 along the whole age range (1 - 13 Gyr). The result¬ 
ing mean dispersion is slightly smaller than typical age er¬ 
rors (A(log(age yr _1 ) = 0.10-0.15) derived from isochrone 
fitting to good-quality cluster CMDs, regarless of binarity, 
multiple populations, differential reddening, rotation effects, 
etc. This result shows that the age-metallicity diagnostic di¬ 
agram returns accurate ages at low expense, since it does 
not require deep photometry nor deal with the known 4- 
parameter degeneracy (age, metallicity, distance, redden¬ 
ing) when matching isochrones to the cluster CMDs. On 
the metallicity arena, the derived iron to hydrogen ratios 
are within A([Fe/H]) = ± 0.15 dex of the identity relation 
(see Fig. 7), no systematic dependence with the metallic¬ 
ity taken from the literature is visible along the considered 
range. Such a dispersion is also comparable to the smallest 
error atainable in deriving metal abundances from photo¬ 
metric data. The most discrepant point in Fig. 7 corresponds 
to NGC6791, which is simply due to a drop in the metallic¬ 
ity sensitivity for metal-rich clusters older than 5 Gyr. 

These results reveal that the diagnostic diagram is able 
to solve the constrains found in the STi index, namely, the 
loss of age sensitivity for ages greater than ~ 6 Gyr and 
the strong dependence on both age and metallicity for ages 
smaller than ~ 5 Gyr. Indeed, as for metal-poor and old 
clusters, the loss of age sensitivity in 5T\ is surpassed when 
the SC - STi index is used as a variable instead of the age as 
proposed by G97 (see top left-hand panel of Fig. 2). Such a 
choice also allows to obtain a metallicity estimate. For ages 
younger than ~ 5 Gyr, the diagnostic diagram shows that 
while STi depends on the age and on the metallicity, SC - 
STi is mainly a metal abundance indicator, so that the latter 
fixes the metallicity level where the former is evaluated. 

In addition to the 5T\ age calibration, the Mr, versus 
(C — Ti) d CMD was also calibrated in terms of metallic¬ 
ity by jCeisler fc Sarajedinij (1999), who demonstrated the 
metallicity sensitivity of the standard giant branch (SGB, 
each giant branch corresponds to an iso-abundance curve) 
applicable to objects with ages > 5 Gyr, any age effects are 
small or negligible for such objects. However, the SGBs were 
defined for [Fe/H] < -0.5 dex using globular clusters older 
than 10 Gyr, so that it is important to examine as closely 
as possible the effect of applying such a calibration to much 
younger clusters. In view of the well-known age-metallicity 
degeneracy, |Bica~et al. (1998| explored this effect empirically 
by comparing SGB-based metallicities for 11 clusters with 
ages between 1 and 3 Gyr to standard values. They found 
a relatively constant offset of ~ 0.4 dex, the SGB metallici¬ 
ties were underestimated due to the effects of age for clusters 
younger than ~ 3 Gyr. Geisler et al. (2003) investigated this 


effect in much more detail by using theoretical isochrones 


computed by Lejeune & Schaerer (2001) for two metallicity 
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levels ([Fe/H] = -1.3 dex and -0.7 dex). They found that 
not only a constant offset of ~ 0.4 dex but an exponential 
correction increasing towards younger ages is necessary. Par¬ 
ticularly, they adopted the theoretical prediction for [Fe/H] 
= -0.7 dex as the correction to be applied to the SGB metal- 
licities as a function of age. Including all error sources, the 
corrected [Fe/H] values are estimated with an uncertainty 
of cr([Fe/H]) = 0.3 dex, although the steepness of the age 
correction for the youngest clusters (< 2 Gyr) results in a 
larger metallicity error and bias the resulting metallicities 
upwards. 

The referred successive improvements seem surpassed 
by the new age-metallicity diagnostic diagram. The latter 
can be used without the knowledge of the cluster distance 
and reddening, which is mandatory in the SGB technique. 
Likewise, the diagnostic diagram allows not only to estimate 
directly more precise metallicities but also ages, simultane¬ 
ously. The new procedure is subject of no correction and 
useful for a wide range of ages (1 - 13 Gyr) and metallicities 
([Fe/H] = -2.0 - +0.5 dex). All these feature make the diag¬ 
nostic diagram a powerful tool for estimating accurate ages 
as well as metallicities. 


5 CONCLUSIONS 

The Washington photometric system has long been used to 
estimate ages and metallicities of clusters, particularly for 
those older than ~ 1 Gyr. Nevertheless, these estimates have 
relied on calibrations (STi and SGB methods) which involve 
clusters with ages and metallicities that need to be updated. 
Likewise, the well-known age-metallicity degeneracy has not 
been properly addressed or even not considered at all. 

We have profusely exploited synthetic T\ versus C — Tj 
CMDs with the aim of improving our knowledge about the 
intrinsic behaviour of the 5T\ index with age and metallic¬ 
ity. The synthetic CMDs were produced through the ASteCA 
suit of functions, taking into account the total cluster mass 
as a function of age in order to have the MSTO and the 
RC similarly populated to those known low mass clusters. 
Photometric errors were also considered, so that the result¬ 
ing cluster CMDs achieved the appearance of the observed 
ones. 

The analysis of the 5T\ index as a function of the age 
for different metallicity levels shows that it varies with age 
and metal content as well. In general, the dependence on age 
weakens for ages greater than ~ 6 Gyr, and results even less 
sensitive to age as the metallicity decreases ([Fe/H] < -1.0 
dex). For ages younger than ~ 5 Gyr 5Ti shows a strong 
correlation with both age and metallicity. As expected, the 
SC index -defined as 5T\ for C the passband - is also a 
combined measurement of age and metallicity. 

We introduce a new age-metallicity diagnostic diagram, 
STi versus SC - 5T \, which has shown the ability of unam¬ 
biguously providing age and metallicity estimates, simulta¬ 
neously, within certain Washington CTj photometric error 
limits. The proposed technique does not require any addi¬ 
tional measurement from other Washington passbands, but 
only the same CT i photometry needed to measure the for¬ 
mer ST\ index. The new procedure allows to derive ages 
and metallicities within a considerable wide range (age: 1 
- 13 Gyr, [Fe/H]: -2.0 - +0.5 dex), and is independent of 


the cluster reddening and distance modulus. It does solve 
the constraints found in the STi index and surpasses the 
performance of the SGB method. 
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Table 1. Estimated age (Gyr) and metallicity (dex) errors* 


[Fe/H]/ age 

1.0 

2.0 

3.2 

4.0 

5.0 

6.3 

-2.0 

0.20/0.15 

0.20/0.20 

0.30/0.20 

0.50/0.25 

0.70/0.20 

0.80/0.25 

-1.5 

0.20/0.15 

0.20/0.15 

0.25/0.20 

0.50/0.25 

0.60/0.25 

0.70/0.25 

-1.0 

0.20/0.15 

0.30/0.15 

0.40/0.25 

0.60/0.25 

0.60/0.25 

0.60/0.25 

-0.5 

0.20/0.15 

0.25/0.15 

0.40/0.25 

0.50/0.20 

0.50/0.25 

0.50/0.25 

0.0 

0.15/0.20 

0.15/0.20 

0.20/0.20 

0.40/0.20 

0.40/0.20 

0.50/0.40 

+0.5 

0.10/0.20 

0.10/0.20 

0.20/0.20 

0.30/0.20 

0.30/0.20 

0.40/0.50 

[Fe/H]/ age 

7.1 

8.0 

9.0 

10.0 

11.2 

12.6 

-2.0 

0.60/0.20 

0.60/0.15 

0.70/0.15 

0.70/0.15 

0.70/0.15 

0.70/0.15 

-1.5 

1.00/0.20 

1.00/0.15 

1.00/0.15 

1.00/0.15 

1.00/0.15 

1.00/0.15 

-1.0 

0.80/0.25 

1.00/0.25 

1.00/0.25 

1.00/0.20 

1.00/0.15 

1.00/0.15 

-0.5 

0.50/0.20 

0.70/0.20 

0.70/0.20 

0.60/0.20 

0.60/0.20 

0.60/0.20 

0.0 

0.50/0.40 

0.50/0.40 

0.50/0.40 

0.50/0.40 

0.50/0.40 

0.50/0.40 

+0.5 

0.80/0.50 

0.80/0.50 

0.80/0.50 

0.80/0.50 

0.80/0.50 

0.80/0.50 


* Errors in age/metallicity were estimated using cr(SC) = a(5T±) = 0.05 mag and a(5C - 8T\) = [(crC) 2 + (crTi) 2 ] 1 / 2 = 0.07 mag. 


Table 2. Cluster parameters taken from the literature used for comparison purposes. 


Cluster 

Age (Gyr) 

Ref. 

[Fe/H] (dex) 

Ref. 

CT\ data 

Age (Gyr) 

[Fe/H] (dex) 






Ref. 

(this 

work) 

47Tuc 

13.1 ± 0.9 

1 

-0.75+ 0.04 

11 

15 

13.0 ± 1.0 

-0.80 ± 0.30 

AM 3 

4.9 ± 1.8 

2 

-0.8 ± 0.4 

2 

16 

4.5 ± 0.7 

-0.75 ± 0.40 

ESO 121-SC03 

9.0 ± 0.7 

3 

-0.93 ± 0.20 

12 

17 

9.0 ± 1.0 

-0.90 ± 0.30 

HW 40 

2.50 ± 0.35 

2 

-0.90 ± 0.15 

2 

16 

3.0 ± 0.7 

-1.00 ± 0.40 

IC 2146 

1.55 ± 0.05 

4 

-0.4 ± 0.2 

4 

18 

1.5 ± 0.3 

-0.60 ± 0.25 

Lindsay 3 

1.2 ± 0.3 

2 

-0.40 ± 0.15 

2 

19 

1.3 ± 0.2 

-0.50 ± 0.15 

Lindsay 113 

4.0 ± 0.7 

2 

-1.24 ± 0.11 

2 

20 

3.5 ± 0.3 

-1.20 ± 0.30 

NGC 339 

6.0 ± 0.5 

5 

-1.08 ± 0.12 

5 

16 

5.2 ± 1.0 

-1.00 ± 0.40 

NGC419 

1.4 ± 0.2 

6 

-0.67 ± 0.12 

13 

21 

1.3 ± 0.3 

-0.75 ± 0.30 

NGC 2682 

4.2 ± 0.2 

7 

0.00 ± 0.06 

23 

15 

4.2 ± 1.0 

-0.10 ± 0.30 

NGC 6791 

7.0 ± 1.0 

8 

0.42 ± 0.05 

23 

15 

7.5 ± 0.8 

0.00 ± 0.50 

SL509 

1.2 ± 0.2 

9 

-0.54 ± 0.09 

9 

17 

1.0 ± 0.2 

-0.40 ± 0.15 

SL862 

1.7 ± 0.2 

9 

-0.47 ± 0.10 

9 

17 

1.7 ± 0.4 

-0.50 ± 0.25 

Trumpler 5 

3.0 ± 1.0 

10 

-0.40+ 0.05 

10 

22 

3.0 ± 0.5 

-0.40 ± 0.15 


Ref.: (l) |Roediger et al.| (|2014[ > ; ( 2) |Dias et al.| (|2014|>; (3)|Mackey, Payne Sz Gilmore| ( |2006[); (4)|Milone et al.| <|2009[); (5)|Glatt et al.| 
( |2011[ ) ; (6)|Glatt et al.|(|2009| >; (7 ) |Balaguer-Nuhez, Galadi-Ennquez Sz Jordi| (|2007[ ); (8 ) | Anthony-Twarog, Twarog &; Mayer|(|2007[);"(9) 

_ |Sharma et al.| (|2010[ >; ( 10)|Donati et al.| (|2015| ); (1 1) |Bragaglia et al.| (|2010|>; ( 12) | Olszewski et al.~ ( |1991|> ; (13) |Da Costa &; 

Hatzidimitriou| ( |1998|>; (14) |Suntzeff et al.||1992[); (15) Geisler Sz Sarajedini| (|1999| ); (16 ) |Piatti| (|2011a); (17)|Bicaet al.|(1998[ ); (18) 


Piatti| (|2011b|); (19)|Piatti et al. (2011); (20) Piatti et al. (2007|); (21)|Piatti ( 201 lc[); (22) |Piatti, Claria & Ahumada (2004); (23) |Heiter| 

et al.|f2014|). 
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Figure 1. Examples of synthetic CMDs produced using the ASteCA suit of functions. 
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Figure 2. Relationship of SC and ST\ indices with the age for different [Fe/H] (dex) values (left-hand panels) and with the metallicity 
for different ages (Gyr) (right-hand panels). The black curve in the upper left-hand panel corresponds to the G97’s calibration. 
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Figure 3. Synthetic age range to cr(age( 397 ) ratio as a function of 8T\. 
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Figure 4. ST\ versus SC - ST\ diagram with iso-age lines and iso-metallicity locii. Metallicity and age labels are given in dex and Gyr, 
respectively. 
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Figure 5. Same as Fig. 2 (upper left-hand panel), where thick line sections correpond to ages/metallicities with uncertainties in boldface 
characters in Table 1. 
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Figure 6. Comparison between ages taken from the literature and estimated using the age-metallicity diagnostic diagram. 
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Figure 7. Comparison between metallicities taken from the literature and estimated using the age-metallicity diagnostic diagram. 




































